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We derive the low-energy effective field theory of the extended t-J model in the regime of light 
doping. The action includes a previously unknown Berry phase term, which is discussed in detail. 
We use this effective field theory to calculate spin-wave excitations in the disordered spin spiral 
state of La2- :E Sr :E Cu04 (the spin-glass phase). We predict an excitation spectrum with two distinct 
branches: The in-plane mode has the usual linear spin-wave dispersion and the out-of-plane mode 
shows non-trivial doping dependent features. We also calculate the intensities for inelastic neutron 
scattering in these modes. 
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I. INTRODUCTION 

The phase diagram of the prototypical cuprate su- 
perconductor La 2 _ 2; Sr 2 .Cu04 (LSCO) shows that the 
magnetic state changes tremendously with Sr doping. 
The three-dimensional antiferromagnetic (AF) Neel or- 
der identified 1 - below 325 K in the parent compound dis- 
appears at doping x w 0.02 and gives way to the so- 
called spin-glass phase which extends up to x w 0.055. 
In both, the Neel and the spin-glass phase, the system 
essentially behaves as an Anderson insulator and exhibits 
only hopping conductivity. Superconductivity then sets 
in for doping x > 0.055, see Ref. @- One of the most 
intriguing properties of LSCO is the static incommen- 
surate magnetic ordering observed at low temperature 
in elastic neutron scattering experiments. This ordering 
manifests itself as a scattering peak shifted with respect 
to the antiferromagnetic position. Very importantly, the 
incommensurate ordering is a generic feature of LSCO. 
According to experiments in the Neel phase, the incom- 
mensurability is almost doping independent and directed 
along the orthorhombic b axis^. In the spin-glass phase, 
the shift is also directed along the b axis, but scales lin- 
early with doping^£. Finally, in the underdoped super- 
conducting region (0.055 < x < 0.12), the shift still scales 
linearly with doping, but it is directed along the crys- 
tal axes of the tetragonal lattice^. Apart from measure- 
ments at x = 0.024 reported in Ref. [5j, all inelastic stud- 
ies on LSCO have so far been performed at sufficiently 
high doping£i2Jii, where the material is conducting and 
the incommensurate structures are thus always aligned 
along the Cu-0 bonds. Similar incommensurate features 
have also been observed in inelastic neutron scattering in 
YBCO, see, e.g., Refs. or Ref. H3 for a 

review. 

The observation of these static incommensurate peaks 
in LSCO caused a renewal of theoretical interest in the 
idea of spin spirals in cuprate a 18 i 19 ' 20 i 21 ' 22 i 23 . While 
the static spiral in the conducting state (x > 0.055) 



is probably an oversimplification, we strongly believe 
that it represents the right picture of insulating LSCO 
(x < 0.055), where mobile holes are trapped in the vicin- 
ity of Sr dopants. Based on a fully controlled solution of 
the t-J model, we first analyzed the low-doping regime 
(x < 0.02) and demonstrated 2 ^ that local static spirals 
are present in the Neel state. In a recent paper—, we 
extended this model to the spin-glass phase of LSCO 
(0.02 < x < 0.055) and found that the ground state is not 
really a spin glass, but a disordered spin spiral state. To 
prevent confusion, we nevertheless comply with the usual 
terminology and refer to this regime as the spin-glass 
phase. The structure factors obtained in our approach 
are in perfect agreement with elastic neutron scattering 
experiments. 



The effective Lagrangian for the weakly doped t-J 
model has been derived quite some time ago 26 ' 27 ' 28 . The 
kinematic structure of the static limit established in 
Ref. can be used without modifications to calculate the 
above mentioned ground state properties of LSC O 24 ' 25 , 
the only difference being the material dependent param- 
eters. However, the derivation of the time-dependent 
terms, necessary to describe excitations, is non-trivial 
and we thus discuss it in detail in the present work. The 
rest of the paper is organized as follows: In Sec. [HI we de- 
rive the general effective field theory describing the low- 
energy sector of the the lightly doped t-J model. This 
effective action applies to the uniformly doped conduct- 
ing state as well as the disordered insulating state. In 
Sec. the effective model is used to calculate the spec- 
trum of spin- wave excitations and the neutron scattering 
cross sections in insulating LSCO. These predictions can 
be directly verified in inelastic neutron scattering exper- 
iments. Our conclusions are presented in Sec. IIVI 
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II. EFFECTIVE LOW-ENERGY ACTION OF 
LIGHTLY DOPED TWO-DIMENSIONAL 
t — t' — t" — J MODEL 

Over a decade ago, the two-dimensional t-J model 
has been suggested to describe the essential low-energy 
physics of high-T c cuprate o 29 ' 30 ! 31 . In its extended ver- 
sion, this model includes additional hopping matrix ele- 
ments t' and t" to next-nearest neighbors. The Hamilto- 
nian of the model is well known, see, e.g., Ref. 0, and 
we do not present it here. The numerical values of the 
parameters of the t-t'-t"-J model for LSCO follow from 
Raman spectroscopy 3 - 2 - and ab-initio calculations 33 . It is 
convenient to measure energies in units of J w 130 meV, 
so that t = 3.1, t' = —0.6 and t" — 0.3. At zero doping 
(no holes), the t-J model is equivalent to the Heisenberg 
model and describes the Mott insulator La2CuC>4 (LCO). 
The removal of a single electron from this Mott insula- 
tor, or in other words the injection of a hole, allows the 
charge carrier to propagate. Single-hole properties of the 
t-J model are well understood 3 ^. The main features are 
a very flat dispersion along the edges of the magnetic 
Brillouin zone (MBZ) with four degenerate half-pockets 
centered at S = (±f ,±f). The quasi-particle residue Z 
at the minimum of the dispersion is Z pa 0.3. In the full- 
pocket description, where two half-pockets are shifted by 
the AF vector Qaf = (tt, tt), the two minima are located 
at S a = (f , f ) and S b = (§, -f ). The system is thus 
somewhat similar to a two- valley semiconductor. 



A. Static limit 

The relevant energy scale for small uniform doping at 
zero temperature is of the order of xJ -C J. In the case of 
nonuniform doping, when holes are trapped near Sr ions 
in hydrogenlike bound states, the relevant energy scale is 
the bin ding energy of the hole, which is about 10 meV, 
see Ref.ull Since in both situations, the energy scales are 
much smaller than J, one can simplify the Hamiltonian 
of the t-J model by integrating out all high-energy fluc- 
tuations. For the hole-doped case, the effective energy 
density reads 



elastic energy of the nonlinear a model (NLSM) and the 
last term is the interaction of the hole with the twist 
of the n- field that favors the formation of local spirals. 
The kinematic structure of this term was first derived in 
Ref. l2Sj, and the coupling constant g i=a Zt sa J was calcu- 
lated in Ref. The second term represents the single- 
particle energy of a hole. The kinetic energy is expanded 
around the center of the corresponding hole pocket. In 
general, the kinetic energy is anisotropic. However, for 
the values of t', t", and t relevant for LSCO, the dis- 
persion is practically isotropic with j3 sa 2 J, see Ref. 
Note that we set the lattice spacing of the square lattice 
equal to unity. Substituting the physical values for the 
lattice spacing, one can check that (3 sa 2 J corresponds 
to an effective mass of about two electron masses. The 
bottom of the hole band for a given flavor is equal to A a 
and we have the freedom to set Ab = 0. For the tetrago- 
nal structure, the hole pockets S a and Sb are degenerate 
and therefore A a = Ab = 0. This degeneracy is lifted 
by the orthorhombic distortion in the low-temperatire 
phase of LSCO. The difference in energy is simply due to 
the slightly different distances between neighboring sites 
along the orthorhombic a and b directions. According to 
a recent calculation 36 , A Q ~ 7 meV. In the spin-glass 
phase of LSCO, holes therefore only occupy the Sb hole 
pocket, which in turn explains the pinning of the incom- 
mensurate diagonal spin structure to the orthorhombic 
b axis, see Ref. Finally, V(r) is the Coulomb poten- 
tial due to the Sr ions that leads to hydrogenlike bound 
states. In the effective energy (p}, we have neglected 
the Dzyaloshinskii-Moriya and the XY anisotropics, see 
Ref. [33, because these anisotropics are small and there- 
fore not important for purpose of the present work. A 
more formal discussion of the effective energy (JlJ can 
be found in Ref. [38|. We note that the pure t-J model 
(£' = t" = 0) is unstable with respect to phase sepa- 
ration and/or short-range charge stripe formation^ 9 -. In 
this case, the long-wavelength approximation considered 
in this section is meaningless. However, the extended 
t-t'-t"-J model at physical values of t' and t" is stable 
in the charge sector— and the long-wavelength action is 
thus well defined. 



£ = ^(vn) 2 + ]T hi 



/3V 2 



A Q + V(r) 



-V2g{i>la-4> a )-[nx (e Q -V)n]} . (1) 

Here n(r) is the staggered component of the copper spins 
and p s sa 0.18 J is the spin stiffness. if> a ( r ) is a fermionic 
spinor operator describing the holes, with an index a = 
a, b (flavor) indicating the location of the hole in momen- 
tum space (either in pocket S a or Sb)- The operator a 
is a psudospin that originates from the existence of two 
sublattices at zero doping and e Q = (l/y2, ±1/V2) is 
a unit vector orthogonal to the face of the MBZ where 
the hole is located. The first term in Eq. fT]) is the usual 



B. Derivation of the time-dependent Lagrangian 

The static limit, i.e., the effective energy given by 
Eq. (P) , is sufficient to determine ground state properties 
of insulating LSCO 24 ' 25 . In order to calculate excitation 
spectra, one has to include the time dependence of the 
staggered field n and the hole operator rp. Naively adding 
the usual time derivatives, we find the Lagrangian 



L = ——n 



Y.{\^«-^m\-£ , (2) 



where \± ~ 0.5/(8 J) is the perpendicular magnetic sus- 
ceptibility and £ is given by Eq. ([!}. However, due to 
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FIG. 1: (Color online), (a) Characteristic ground state configuration of a particular realization at doping x = 0.05. The 
impurity pseudospins ^Sia^i (circles) are oriented along the z axis. Full (open) circles correspond to values — i (+■§)• Small 
arrows represent the n-field. The system forms domains stretched along the a direction, in which all the pseudospins are aligned 
in parallel, (b) Uniform spiral state (I12|l with average pitch (|13|) . Apart from the absence of disorder, a uniform spiral is very 
similar to the true ground state and thus a good starting point to describe excitations. The corresponding structure factors 
are shown in Fig. [2] 



the interaction of the holes with the n-field [last term 
in Eq. ([T])], there is an additional non-trivial term. Let 
us find its kinematic structure by following the original 
calculation of the static Umrfc2£ and extend it to the time- 
dependent case: The hopping term between nearest sites 
i and j in the t-J model is given by 

H t = -tJ2 4 s c ]S + H.c. , 

where c\ s creates an electron on site i with spin projec- 
tion s = ±1/2 onto a globally fixed z axis. For the Neel 
state with a single injected hole \h^ a ) on the £ =|, | sub- 
lattice, with flavor index a = a, 6, the expectation value 
(hia\ H t \hi a ) vanishes because spins on nearest sites are 
exactly opposite. Let us now assume that the spin back- 
ground is deformed in such a way that the staggered mag- 
netization is given by ft = (sin#cos<y9,sin#sin</?,cos#) 
where and ip smoothly depend on r. Because we assume 
the deformation of the background to be smooth, the two 
sublattices j and J, are still well defined, but the arrows 
no longer represent the projection of the spin along the 
z-axis. After the unitary transformation of the operators 
Ci S into the local reference frame with quantization axis 
z directed along ft, we find 

H t -V2t (h\ a (r)[6' + Vsin0]ft ia (r)+H.c.) , (3) 

where F 1 = (e a ■ V)F. Due to high-energy quan- 
tum fluctuations, t is replaced by t — ► g — Zt, see 
Ref. HH. The Hamiltonian §5§ leads to the following 
Schrodinger equation for the two-component wave func- 



tion v a = (/i Tq , h la ) 

i v a = -V2g(f ■ a)v a , (4) 

where / = (9', tp'sin.6, 0) and a are the Pauli matrices. 
Since 

f = {6') 2 + (ip') 2 sin 2 $=[ftx (e Q ■ V)n] 2 , 

one can perform a unitary transformation U , ip a = Uv a , 
such that 

UHf- a)U = [ft x (e Q ■ V)n] • a . 

The Schrodinger equation (T4|) for a single hole is thus 
transformed into 

i ip a = —V2g [ft x (e Q ■ V)n] • aip a — U^Uipa . (5) 

If the n-field is static, then U = and the Schrodinger 
equation ([5]) is equivalent to the Euler-Lagrange equa- 
tion that follows from the effective action ([2]). Taking 
into account the U contribution in {5| thus leads to the 
introduction of an additional term in the effective ac- 
tion. Its kinematic structure can be determined due to 
the following reasons: 1) The term originates from Wll 
and thus contains only one time derivative ft. 2) Because 
it acts on a pseudospinor, it contains a. 3) It must be 
a scalar. 4) It is independent of the coupling constant 
g. The only structure that satisfies these conditions is 
i/j^ffi/ja ■[Hx ft]. A calculation shows that it appears with 
a universal coefficient —1/2 and therefore leads to the 
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Lagrangian 



L = ——n 



rv *• 



I ■ [n x n}\ - £ . (6) 



The additional a- [fix ft] term, closely related to the Berry 
phase, has been obtained before for Kondo-lattice-like 
model o 39 ' 40 . However, in these models the term appears 
with a non-universal coefficient that depends on the pa- 
rameters of the microscopic description. It is important 
to note that the effective Lagrangian ((6]) is not restricted 
to the situation of localized holes, but is also valid in the 
uniformly doped case when V(r) in Eq. (JTJ is zero. 



C. Larmor's theorem 

The validity of the prefactor — | in the Berry phase 
term can be verified in a situation where the behavior 
of the system is well known: According to Larmor's the- 
orem, spins in a uniform magnetic field B process with 
frequency lj = B. In the undoped case, this property 
follows directly from the Euler-Lagrange equations. Ac- 
cording to Refs. I4lll42l . the Lagrangian of the NLSM in 
an external magnetic field is given by 



— ( ft — [ft x B 
2 1 



|[Vn] 2 



Note that we set g s pB — * 1, gs being the gyromagnetic 
factor and fiB the Bohr magneton. For a uniform field 
Vn = 0, the Euler-Lagrange equation ft = [ft x B] de- 
scribes the precession of the staggered field around B, in 
accordance with Larmor's theorem. 

In a doped system, the situation is similar. Let us 
for simplicity consider a single hole trapped by the Sr 
potential U(r) and omit the flavor index a. The general- 
ization to multiple dopants is straightforward. Using the 
pseudospin density 



we perform the variation of Eq. 
and find 



with respect to n(r) 



A? 



+ ^ f [ft x f] + 2[n' x |f) + An = . (7) 



Here F' = (e Q -V)F and A = A(r) is a Lagrange multiplier 
taking into account the constraint n 2 = 1. It has been 
shown in Ref. |2l| that n(r) that results from Eq. ([7]) lies 
in a plane to which the pseudospin is perpendicular. In 
other words, £ is parallel to ft x ft' . 

If we place the doped system in a uniform magnetic 
field, we have to replace 



[ft — [ft x B]j 



in the Lagrangian ([6]), and we also have to take into ac- 
count the interaction energy of the pseudospin with the 
magnetic field, Sr = — % (£ ■ n){B ■ ft), that has been de- 
rived in Ref. [24j. In the presence of a uniform magnetic 
field, the effective Lagrangian thus reads 



C B = ^ ( ft - [ft x B 



- y (Vn) 



+ o(V ,t V' - flip) - ^ 



/3V 2 



V(r) 



-I*- [fix ft] + V2g£ ■ [ft x ft'} + -(I*- ft)(B ■ ft) 



(8) 



Here we have omitted the long derivative p — > p — |A, 
describing the interaction of the magnetic field with the 
electric charge, because we are only interested in the spin 
dynamics. Larmor's theorem implies that ft — ft x B 
still satisfies the equation of motion. Substituting this 
solution into Eq. ©, we find 



£u--f [Vn] 2 - 



- — ( 1p ^ 1p — lj) ^ "0 ) — "0 1 



/3V 2 



V(r) 



1> 



v / 2<3 , £*- [n x ft'] + -(£*• B) . (9) 



Since according to Eq. (O, ^ is parallel to ft x ft', the 
equation of motion that follows from the Lagrangian ([9]) 
is 



(10) 



i.e., £ also precesses with frequency lj = B. Eq. (J7J) there- 
fore remains valid in the proper reference frame and we 
conclude that the Berry phase term is crucially important 
for the fulfillment of Larmor's theorem. 



III. SPIN- WAVE EXCITATIONS IN THE 
SPIN-GLASS PHASE OF LSCO 

In the Neel and the spin-glass phase, holes are local- 
ized near Sr ions in hydrogenlike bound states 2 ^. In what 
follows, we refer to these bound states as impurities. 
The binding energy is about 10 meV and the descrip- 
tion is therefore valid at temperatures well below 100 K. 
A direct experimental indication of this picture is the 
variable-range hop ping conductivity observed in insulat- 
ing LSCO, see Ref. SI or Ref. Q for a review. 



A. Uniform spiral approximation for the ground 
state 

In the case of hole localization, the hole wave function 
ipi(r) decays exponentially away from the Sr ion located 
at position 
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FIG. 2: (Color online), (a) Neutron scattering probability Sq for doping x = 0.04. Full circles correspond to experimental 
observations taken from Fig. 4 in Ref. [(J with normalized intensities. The solid line represents our simulation^, containing 
no fitting parameters. The dashed line indicates the positions of the delta function peaks obtained from a uniform spiral, (b) 
Incommensurability S (half the distance between the two peaks measured in reciprocal lattice units of the tetragonal lattice), as 
a function of doping. Accurate Monte Carlo calculations (dots) are in agreement with experimental data (squares) taken from 
Refs. OHHm, see Fig. 6 of Ref. H. The uniform spiral (dashed line) captures the essential properties of the spin-glass ground 
state. 



The energy ([T]) is then given by 



£ = — [Vn] - 



V2g^2 P ( r - R i) \n x ( e & ' v ) A > 

i 

(ii) 

with p (r) = §t 2 e 2kt and the inverse localization length 
k ~ 0.4 that follows form hopping conductivity^. We re- 
mind the reader that we set the lattice spacing of the 
tetragonal structure equal to unity. In absolute units, 
the hole localization length is thus approximately 10 A. 
Note that the pseudospinor of a given impurity i is in- 
dependent of r. We do not include the binding energy in 
(fTTjl because it just shifts the energy by a constant. The 
ground state is found by minimizing the energy . Af- 
ter integrating out the n-field, this minimization proce- 
dure reduces to finding the minimum of a system of ran- 
domly distributed interacting dipoles & = ^i<j^i. This 
task has been accomplished in our previous work>2£ using 
Monte Carlo simulations. A characteristic ground state 
configuration of the dipoles and the n-field obtained for 
a particular realization of Sr positions at average doping 
x = 0.05 is shown in Fig. [IJa). We would like to empha- 
size that this state is not really a spin-glass, but a disor- 
dered spin spiral stated. The spiral structure manifests 
itself as incommensurate peaks in the structure factor 



V (r,)n Q (pj) 



Our claim about the disordered spiral state is supported 
by experimental observations: While the incommensu- 
rate structure is already observed below 30 — 40 K, the ir- 
reversible glassy behavior only sets in below T ~ 5— 6 K 44 
and is in our opinion due to the inter-layer interaction 
that leads to a freezing of incompatible spiral configu- 
rations in different planes. The profile of the neutron 
scattering cross section calculated in Ref. is shown 



in Fig. Ufa), together with experimental data at dop- 
ing x — 0.04. Our results are in good agreement with 
experiments, especially given the fact that they contain 
no fitting parameters. Even though the values for the 
important parameters k and g are fixed, there is some 
variability related to their accuracy. In the present work, 



as well as our two recent paper 



,24,25 



we use k 



0.4 



and g « 1. However, we believe that k = 0.3 or g = 0.8 
are well possible. The spiral has an average pitch propor- 
tional to doping, as can be seen in the doping dependence 
of the incommensurability S, defined as half the distance 
between the two peaks, shown in Fig.[^b). 

The observed broadening of the incommensurate peaks 
is on the one hand due to intrinsic disorder caused by the 
random distribution of the Sr ions and on the other hand 
due to spin vortices (topological defects) that determine 
the size of the domains with a given chirality of the spiral, 
see Ref. [25l These broadening mechanisms can only be 
taken into account in numerical studies. However, if we 
neglect these effects, we can find the ground state in a 
much simpler mean-field approach: Since the interaction 
term in Eq. pip is maximal when all the pseudospins are 
orthogonal to plane of the spiral, we can take 



& = *ia% = (0, 0, 1) , 

n = (cos Q ■ r, sin Q ■ r, 0) . 



(12) 



After averaging over the impurity positions, we find the 
energy £ = p s /2Q 2 — \/2gx(ei, ■ Q). The spiral pitch that 
minimizes this energy is given by 



Q = e h , 



(13) 



and agrees well with accurate Monte Carlo simulations 
and experimental data, see Fig. [U Fig. [2] and Ref. HE]. 
Hence, apart from the absence of disorder, a uniform 
spiral (fT2")) with average pitch (fT3"|) is very similar to the 
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k b /(2jt) 



FIG. 3: Spectrum of in-plane spin waves excited in neutron 
scattering at doping x = 0.05. The momentum is directed 
along the orthorhombic b direction and the incommensurable 
vector (|13[) is Q « 0.39. There are two linear branches that 
start at kb = ±Q, that are broadened due to disorder and 
ultimately disappear at |k ± Q| ~ \fx. For larger momenta, 
spin waves reappear with usual dispersion of the Neel anti- 
ferromagnet. The frequency at the intersection k(, = scales 
linerly with doping. 

true ground state and is thus a good starting point to 
describe excitations. 



B. In-plane spin-wave excitations 

An in-plane excitation is described by a small deviation 
(p = ip(t,r) from the uniform spiral ground state (|12[) . 
Substituting the expression for the n-field 

n = (cos(Q • r + <p), sin(Q ■ r + </?), 0) . (14) 

into the effective Lagrangian j6]), we find 

£ = ^ 2 -f (Q 2 + 2Q^' + (V^) 2 ) 

+ ]Tp(r - Rj) <.V2g&(Q + if') - ~&y> 

+ ~ , (15) 

where £j = ^ja^i. It should be clearly understood that 
Eq. (fT4"]) and the Lagrangian (JTSJ) are only valid in the 
long- wavelength limit, q < y/x, where q is the typical mo- 
mentum of the y-field and \Jx is the inverse average dis- 
tance between impurities. In this limit, one can average 
over impurity positions and replace J2i P( r ~ R*) ~~ * x > 
which leads to the cancellation of terms proportional to 
(p' in the Lagrangian (fTS")) because of Eq. (fT5|) . From 
the Euler-Lagrange equation for 'J, we find £ = 0. The 
equation of motion for ip therefore reads 

Xx<p = PsV 2 ^ . 



The spectrum of in-plane excitations is thus exactly the 
same as the spin- wave spectrum in undoped LCO, u) — 
cq, with the spin- wave velocity c = y/p s /x± ~ 1-66. The 
corresponding Green's function of the ip-Held is given by 

r - xl 1 

ui z — c z q z + iv 

where the subscript in stays for in-plane. 

Let us also calculate the neutron scattering cross sec- 
tion, in order to make specific predictions for future ex- 
periments. The Hamiltonian describing the interaction 
of the neutron spin S N with the n-field reads 

H N oc S N ■ ft = S? n z + i (S*V + S N n + ) . (16) 

After the substitution of the in-plane excitation (fT4")l . the 
above Hamiltonian reads 

2 + 2 ~ 

-> ~e ikr {S^V iQr (l - tip) + S w e lQ ' r (l + tip)} , 

where k is the momentum transfer and Q the momen- 
tum shift due to the spiral ground state. The scattering 
probability for unpolarized neutrons is given by 

Ji„(w,k) oc -— Im [G m (w,k- Q) + G in (uj,k + Q)] 

= -J_[ ( y(a;- c |k-Q|) + 5(w-c|k + Q|)] . (17) 
16x_i_w 

Fig. [3] shows the spectrum of excited spin waves for a mo- 
mentum transfer k = (0, k^) chosen along the orthorhom- 
bic b direction at doping x — 0.05. Very importantly, the 
above calculation is valid at q = |k± Q| < y/x and in the 
absence of disorder. In experiments, we expect a broad- 
ening similar to that obtained in elastic neutron scat- 
tering, see Fig. [5J For momenta q ~ \fx (shaded area 
in Fig. [3]), the peaks disappear because of very strong 
broadening due to scattering of the spin- wave off random 
impurities. However, we predict that at even higher mo- 
menta, the peaks reappear as usual spin waves, to — ck, 
with some broadening, reflecting the nearly perfect anti- 
ferromagnetic alignment of nearest neighbor spins. 

C. Out-of-plane spin-wave excitation 

An out-of-plane excitation can be represented as a 
small deviation n z (r, t) from the uniform spiral ground 
state (fT2")) . Let us write the n-field as 

ft = ( \J 1 — n 2 z cos Q • r, \J\ — n\ sin Q • r, n z ) , 

and substitute this expression into the effective La- 
grangian (J6J. Neglecting cubic and higher order terms 
in n zi we decompose the Lagrangian into two parts 
C = C + A n t, with 
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t 




Q 



FIG. 4: Polarization operator describing the scattering of a 
spin- wave (dashed line) off a pseudospin (solid line). The 
corresponding pesudospin-flip vertices (circles) are given by 
Eq. ®. 



Co = fn^-| {-Q 2 + (Vn z ) 2 + Q%t} + £ 



V2g(Qn z - in' z ) + -n z 



JQr 



V2g(Qn z + in' z ) - -n z 



(18) 



Here A = 2\f2gQ = 4g 2 x/p s is the energy required to flip 
a pseudospin and F' = (e^ • V)F. From the diagonal part 
of the Lagrangian Co, we find the unperturbed spin- wave 
Green's function 



LU 2 — W 2 



(19) 



where cj 2 = c 2 (q 2 +Q 2 ). Note that this is different from 
the usual spin-wave dispersion. In order to calculate the 
corrections to the bare propagator due to Ant, we in- 
troduce the second quantization representation for the 
n-field 



E 

q 



1 



e <u,t-tq.r m t 



with the spin-wave creation and annihilation operators 
ml and m q . The off-diagonal part of the Lagrangian 
Ant (O gives rise to the pseudospin-fiip vertices 

1 



[ I 



IT 



X-L L 2 



1 



2 (Q- q) .R !F(Q _ q) y /2g(Q + qb ) + 



(20) 



In contrast to the in-plane excitations, the Berry phase 
term derived in Sec. Ill Bl is important for the out-of-plane 
spectrum because it leads to the additional terms ±w/2 
in the above expression. The impurity form factor -F(p) 
is given by 



F(p) = J d 2 r p{vyv* = ^-j 



= ^- I d 2 re- 2Kr e tp r 



[1+ P 2 /(Ak 2 )] 



3/2 



8 k 2 



(21) 



where for the last expression, we have used the small 
momentum limit p <C K. 

The polarizability of a given impurity is described 
by the diagram shown in Fig. [4] Using the bare 
Green's function (fT9"|) together with the pesudospin-flip 
vertices (|20|) . the polarization operator is given by the 
polarizability of a single impurity multiplied by the con- 
centration x 



n\ 



ir 



IT I 



A -uj-A 



— {F 2 (Q + q)H(-u>, -q) + F 2 (Q - q)H(u, q)} , 
where 



ff(u,q) 



(Q + q b + Quj/Af 



LU + A 

The polarization operator is invariant under a simulta- 
neous change of the signs of q and lo, i.e., P(u>, q) = 
P(— lu,— q). However, it is neither an even function of 
the frequency lu nor the momentum q. In order to find 
the spectrum for small momenta, we substitute the ap- 
proximate expansion for the form factor (|2ip and neglect 
quartic and higher order terms in q. In this limit, the 
small asymmetry disappears and we find 



P(w,q) 



1-LU 2 /A 2 



An 2 



[(Q 2 -<i) 2 + cil{Q 2 + q l)} 



(22) 



Including single-loop corrections, the Green's function of 
the n, -field reads 



Go 



X± 



j 2 -u;l-P(u;, q ) ' 



8 



300 



200 



100 - 




t/(2it) 



FIG. 5: Spectrum of out-of-plane spin-wave excitations at 
doping x — 0.05. The Momentum is directed along the or- 
thorhombic b direction. The inverse radius of the impurity 
wave-function is k — 0.4 and the incommensurable vector (|13[) 
is Q ss 0.39. Both, the lower and upper branch are getting 
broad and ultimately disappear due to disorder at |k±Q| ~ ft. 
The upper branch has a very small spectral weight and it thus 
difficult to observe in neutron scattering. At larger momenta, 
the spin waves reappear with usual dispersion of the Neel anti- 
ferromagnet. The frequency at the top of the dome at k(, = 
scales quadratically with doping. 



where the subscript out stays for out-plane. The excita- 
tion spectrum is given by poles of the Green's function. 
Using the symmetric expression for the polarization op- 
erator 11221) . we find two branches 



\ 



<& 
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Tk 1 



[(Q 2 -q 2 b ) 2 + q 2 a (Q 2 + q 2 b : 
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The lower branch vanishes at q = ±Q, in accordance 
with the Goldstone theorem, while the upper branch has 
a gap A. The dispersion shown in Fig. [5] is obtained 
from the exact expression of the polarization operator. 
The small asymmetry, clearly visible in this numerical 
calculation, especially in the gapped branch, gives an 
additional broadening after averaging over domains of 
different chirality. Our results are valid for momenta 
q = |k ± Q| < -s/x and in the absence of disorder. In 
this range of momenta, disorder is expected to lead to 
some broadening of the neutron scattering peaks, as ob- 
served in elastic scattering, see Fig. [5J This broadening 
is getting stronger for momenta corresponding to the in- 
verse separation between impurities, where spin waves 
are scattered off random impurities. In this region, the 
peaks are very broad and thus difficult to detect. How- 



ever, we expect a distinct signal to reappear at higher mo- 
menta, reflecting the usual spin-wave dispersion ui = ck 
This reappearance is due to the nearly perfect antiferro- 
magnetic alignment of nearest neighbor spins. The neu- 
tron scattering probability defined by the interaction (fl~6l 
reads 



I out (u>,q) oc -— Im G „ t (cj,q) 

47T 



(23) 



1 



[Zq6(w - Oi, q ) + (i - z q )<j(w - n a , q )] , 



with the quasi-particle residue 



1 



1 + 



Since the upper branch has a practically vanishing spec- 
tral weight, it is difficult to observe experimentally. Com- 
pared to the in-plane spin waves excited in neutron scat- 
tering, the out-of-plane spectrum is very similar, espe- 
cially around momenta q ► Q where the energies go to 
zero. The difference can be best seen around momenta 
q sa 0. However, one should remember that there is con- 
siderable broadening of the peaks due to disorder in this 
region. Interestingly, the ratio between the predicted in- 
tensities for the low-frequency in- and out-of-plane spin 
waves is I out : I in — 2, see Eqs. (fTT)) and (|2"3"j) . Such a 
difference should be detectable in experiments. 



IV. CONCLUSION 

To conclude, we have derived the low-energy effective 
field theory of the t-J model in the limit of small dop- 
ing. Based on this description, we have calculated the 
spin-wave excitations in the disordered spin spiral state 
of I^-^Sr^CuO^ (the spin-glass phase), which has a 
coplanar spin structure. For the in-plane spectrum, we 
have found the usual linear spin-wave dispersion, while 
the out-of-plane modes have been shown to exhibit non- 
trivial doping dependent features. 
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